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ABSTRACT 

The average rate of emission of ionizing radiation per unit volume (or emissiv¬ 
ity) in the universe can be calculated as the ratio of the intensity of the ionizing 
background to the mean free path of ionizing photons. The intensity of the back¬ 
ground is measured from the mean transmitted flux of the Lya forest, and the 
mean free path is measured from the abundance of Lyman limit systems, which 
has been observed so far up to z = 4. This yields an emissivity that is not larger 
than 7 ionizing photons per Hubble time for each atom in the universe a.t z = 4, 
which may reasonably arise from QSOs and star-forming galaxies. In order for 
the reionization to end hj z = 6, and assuming that the clumping factor of ion¬ 
ized gas during the reionization epoch is close to unity, this ionizing comoving 
emissivity cannot decline from 2 ; = 4 up to 2 ; ~ 9 by more than a factor 1.5. 

If the clumping factor were much larger than unity, then the emissivity would 
need to rapidly increase with redshift. Unless the ionizing emissivity increases 
substantially from z = 4 to z ~ 10 — 20, the Thomson optical depth to the CMB 
must be in the range 0.045 < r < 0.09. 

Subject headings: cosmology: theory - diffuse radiation - intergalactic medium 
- galaxies: formation 


1. Introduction 

The highest redshift quasar known at the present time, a.t z = 6.28, is the hrst one to 
show a complete Gunn-Peterson trough (Becker et ah 2001; Djorgovski et ah 2001). This 
suggests that the reionization of hydrogen was probably completed near the epoch z = 6. 
This conclusion is not yet completely certain because an ionized medium can produce a 
complete Gunn-Peterson trough if the neutral fraction is everywhere high enough (owing to 
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a low intensity of the ionizing background). However, the abrupt change with redshift near 
z = 6 of the fraction of transmitted flux in the Lya region of the spectrum (Becker et ah 

2001) implies a rapid increase of the mean free path of ionizing photons with time (Fan et ah 

2002 ) , at least in the region of the line of sight to the = 6.28 QSO, just as expected when 
the last low-density regions of the intergalactic medium are ionized (e.g., Gnedin 2000). 

Ending reionization hy z = 6 requires that at least one ionizing photon has been emitted 
to the intergalactic medium for each baryon in the universe, in addition to the amount 
needed to compensate for any recombinations that have taken place. This implies that the 
emissivity of ionizing photons (dehned as the mean number of photons emitted per unit time 
and volume) should be greater than the mean number density of baryons divided by the age 
of the universe at z ^ 6, unless the emissivity had declined rapidly with time at z > 6. 


It is possible to infer the emissivity at lower redshift, by using the intensity of the 
ionizing background measured from the mean transmitted flux in the Lya forest, and the 
mean free path of ionizing photons determined from the observed abundance of Lyman limit 
systems. The first of these two quantities, the intensity of the ionizing background, can 
be measured in the context of the present theoretical understanding of the Lya forest as 
arising from a photoionized intergalactic medium that follows the gravitational evolution of 
primordial fluctuations (see the review by Rauch 1998 and references therein). The observed 
mean transmitted flux can then be used to infer the Lya scattering optical depth at the 
mean hydrogen density and in Hubble expansion, r^: 
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Here, x^t is the hydrogen neutral fraction at mean density, T 4 is the gas temperature in 
units of 10"^ K, and r_i 2 is the photoionization rate due to the cosmic background in units 
of 10“^^ s“h When the value of Qi,h^ was still highly uncertain, the value of inferred from 
observations of the mean transmitted flux, ~ 1.0 ± 0.2 at z = 3 (see Rauch et al. 1997, 
and Table 2 in McDonald et al. 2000) implied a lower limit to the value of given by the 
lower limit on r_i 2 that was inferred from the observed QSOs. Now, has been more 
accurately measured from the primordial deuterium abundance (O’Meara et al. 2001) and 
the CMB spectrum of fluctuations (de Bernardis et al. 2002; Pryke et al. 2002; Spergel et 
al. 2003), and we can use the value of Tu(z) determined from the Lya forest to constrain the 
intensity of the ionizing background and its evolution with redshift. 


The mean free path of ionizing photons is related to the abundance of Lyman limit 
systems, which have column densities Nffi > 1.6 x lO^^cm”^. The reason is that the column 
density distribution of absorption systems is approximately of the form /(iV^^z) dNjjj oc 
dNni (e.g., Petitjean et al. 1993), so the average rate of absorption is dominated by 
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high column densities, up to the column density at which the absorption systems become 
optically thick. The abundance of Lyman limit systems has been reasonably well measured 
up to redshift z = 4. In this paper, we combine the measurements of the intensity and the 
mean free path of ionizing photons to infer their rate of emission at z = 4 (§2). We then 
compare this emissivity to the observations of known sources at this redshift (§3). Then, in 
§4, we hnd the consequences for the evolution of the emissivity at higher redshift to satisfy 
the requirement of reionizing the IGM hj z = 6. Finally, in §5 the consequences for the 
Thomson optical depth of the Cosmic Microwave Background photons are discussed. We 
will use the flat cosmological model with Qm = 0.3, h = 0.65, and = 0.022. 


2. The ionizing emissivity at z = 4 

In this section we evaluate the emissivity of ionizing radiation at redshift z = 4. In 
terms of the photoionization rate, T, and the mean free path of the ionizing photons. A*, the 
photon emissivity is given by e* = r/(Ajd'), where ct is a frequency-averaged cross section. 
The result we will hnd in this paper is that the emissivity at z = 4 has a rather small 
value, implying that it must not have declined much with redshift, up to a substantially 
earlier epoch at 2 : > 6, in order to satisfy the requirement that reionization was completed 
by ~ 6. With that in mind, we will evaluate the maximum allowed value of T and of at 
z = 4, which is the highest redshift at which the mean free path A^ has been observationally 
determined. 


2.1. The photoionization rate 

The most recent measurements of the mean transmitted flux through the Lya forest, F, 
were reported by Bernard! et ah (2002), from a large sample of QSOs of the Sloan Digital Sky 
Survey. From their Fig. 4, the effective optical depth [dehned as F = exp(—Te//)] at 2 : = 4 is 
Teff = 1.02 ±0.02, where we have estimated the error from the errorbars and scatter of their 
measured points at various redshift bins. Hence, F < 0.37 at z = 4. Note that previous 
values of F from a small sample of high resolution spectra (Schaye et ah 2000, McDonald 
et al. 2000) were higher, but this was probably because their method underestimated the 
QSO continuum flux, which becomes increasingly difficult to set by htting regions between 
absorption lines as the redshift increases. 

We now use the results of McDonald & Miralda-Escude (2001, hereafter MMOl) to infer 
the photoinization rate F_i 2 . According to their Figure 2, the value F < 0.37 implies a 
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photoionization rate of r_i 2 < 0.40 at z = 4.5. This value needs to be corrected from the 
model used in MMOl to the model we use here, and to z = 4, where we will choose all the 
parameters to maximize r_i 2 , and therefore infer an upper limit. First of all, the relationship 
between F and in equation ( 1 ) depends on the amplitude of the mass power spectrum at 
the Jeans scale of the Lyo; forest: a lower amplitude of fluctuations implies that the gas is 
less concentrated in collapsed regions that correspond to saturated absorption lines, and a 
greater fraction of the gas is left in low-density regions, where it contributes more efficiently 
to decrease the mean transmitted flux. The measured amplitude of the power spectrum 
of the Lya forest (see Croft et al. 2002 and references therein) is 0.86 times lower than 
the amplitude in the model that was used by MMOl (as shown in McDonald et al. 2000), 
which had the parameters ag = 0.79 and n = 0.95 for the power spectrum. In other words, 
the model used by MMOl would match the Lya forest power spectrum amplitude found by 
Croft et al. if its normalization had been changed to Ug = 0.86 ■ 0.79 = 0.68, leaving all other 
parameters unchanged. Since the amplitude of mass fluctuations grows proportionally to the 
scale factor a (the effects of the vacuum energy and the radiation are negligible at z = 4), the 
model with normalization ag = 0.68 has the same relation of F and at 2 : = 4 as the model 
with normalization ag = 0.79 at z = {1 + 4)/0.86 — 1 = 4.8. It was also found in MMOl that 
the value of Tu inferred from a hxed F depends on the amplitude approximately as (this 
can in fact be checked by comparing the values of r_i 2 at z = 4.5 and = 5.2 for hxed F in 
Fig. 2 of MMOl). Because the limit F_i 2 < 0.4 was derived at z = 4.5, and F_i 2 oc the 
power spectrum amplitude correction is a factor (5.8/5.5)^. We note that here we have taken 
the central measured value of Croft et al. of the power spectrum normalization, instead of 
a lower value within their errorbar (which would further raise the upper limit to F_i 2 by 
about 20%), because other measurements of the power spectrum amplitude favor a higher 
value instead (see Seljak 2001; van Waerbeke et al. 2001; Croft et al. 2002; Jarvis et al. 2002; 
Spergel et al. 2003). 

Second, the inferred F_i 2 depends on the gas temperature. MMOl assumed T 4 = 2, 
but measurements of the gas temperature from the Lyo; forest spectra have given values in 
the range 1.2 < T 4 < 2 (Theuns et al. 2002 and references therein). Although the optical 
depth of a uniform medium is proportional to (eq. 1 ) because of the temperature 

dependence of the recombination coefficient, the relation between F and Tu is also affected by 
the temperature because thermal broadening spreads some of the absorption from saturated 
spectral regions to other regions with less absorption. We use the result found in McDonald 
et al. (2001, their §7.4) that thermal broadening causes a change in Tu of 27% of that 
due to the recombination coefficient, in the opposite direction, therefore giving an effective 
dependence F_i 2 oc If we assume the lowest value of the temperature among the 

reported measurements, T 4 = 1 . 2 , then the temperature correction for the upper limit of 
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r_i2 is {2/l.2f-\ 

Finally, we correct for all other parameters used in MMOl = 0.02, h = 0.65, 

= 0.4, = 4.5) to the ones used here = 0.022, h = 0.65, = 0.3, z = A), which 

we choose within the current uncertainties to maximize the inferred r_i 2 . The result of all 
these corrections (see eq. 1) for the upper limit to r_i 2 at 2 : = 4 is 
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2.2. The mean free path 

Storrie-Lombardi et ah (1994) found that the rate of Lyman limit systems (dehned as 
absorption systems with column density Nhi > 1.6 x 10^^ cm“^) at 2 : = 4 is dNn/dz = 3.3 ± 
0.6. In other words, the mean spacing between Lyman limit systems is Xll = cH~^{z)/{l + 
z)/{?).?) ± 0.6) = (0.06 ± Q21)cH~^{z) at 2 : = 4. To hnd the relationship of Xu to the 
mean free path of ionizing photons, we assume that the column density distribution of the 
absorbers is /{Nhi) dNni oc dNni (Petitjean et ah 1993). The absorption probability 
per unit length, Aq for a photon at the hydrogen ionization edge frequency u = uhi, is 

1 _ drr-i-5 (1 - e-") _ 

Ao /“drr-i-5Ai^ Xll' 

where r = Nhi/{1-Q x lO^'^cm”^). At a higher frequency u, and approximating the pho- 
toionization cross section as proportional to z/“^, the mean free path is increased to X,y ~ 
Xoliz/ izhi)^'^- 

We note here that our mean free path in (3) is 1.5 times longer than that adopted in 
Madau, Haardt, & Rees (1999), and in Schirber & Bullock (2002). For some reason, their 
opacity corresponds to a Lyman limit system abundance of 5 per unit redshift, larger than 
reported by Storrie-Lombardi et al. (1994). 


2.3. The emissivity 


The photoionization rate is related to the emissivity the cross section and the 
mean free path as a function of frequency by 
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Here, we use the approximation that the mean free path is much less than the horizon 
distance (the exact relation between the emissivity and background intensity involves an 
integration over time with redshift effects, a full knowledge of the spectrum, and reemission; 
see Miralda-Escude & Ostriker 1990; Haardt & Madau 1996; Fardal, Giroux, & Shull 1998. 
We treat reemission approximately in §2.4). Assuming that oc oc and 

(7y oc z/“^, we have T = eoAoCTo/(1.5+ ?]), where ao = 6.3 x 10“^® cm^, and eo is the emissivity 
per unit frequency at the ionization edge. The emissivity integrated over frequency is Cj = 
J due„ = eo/rj. The emissivity can be conveniently expressed as the number of ionizing 
photons emitted per Hubble time and for each atom in the universe: 

= 4) ^ v/i(r/ + 1.5) T 

Ha Uat (0.06 ± 0.01) rjcaoUat (1 + z)^ 

where Uat is the comoving number density of atoms (adding hydrogen and helium), Uat = 
X^bPcrit/'m-Hi^ — 3y/4) = 2.03 x 10“^cm“^, and Ha is the Hubble constant dX z = A. Using 
r_i 2 < 0.53, taking the lower limit for the mean free path (to obtain an upper limit for the 
emissivity), and assuming r] = 1.5, the hnal result is 
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In other words, there are fewer than eight ionizing photons being emitted per Hubble time 
and for each atom at z = 4. 


2.4. Effect of reemission from the Lyo absorbers 

Some of the photons that contribute to the photoionization rate in the intergalactic 
medium arise from reemission by the Lyo absorption systems (see Haardt & Madau 1996). 
For every hydrogen atom that is photoionized by absorbing a photon, and assuming ionization 
equilibrium, a recombination will take place, and about 40% of the recombinations are direct 
to the ground state, producing another ionizing photon. However, the majority of the ionizing 
photons that are not very close to the hydrogen ionization edge are absorbed deep within 
Lyman limit systems, and the photons produced by recombinations in these regions are 
close to the ionization edge [with a typical frequency separation from the Lyman limit of 
Au/ V ~ kT /(13.6 eV) ~ 0.1] and not likely to escape from the Lyman limit system. Because 
of this, reemission by hydrogen recombination increases the background intensity by only 
~ 10% (as is found, for example, by using the approximation that the reemitted photons 
escape only when produced in an absorber of column density N^i < 2 /(To, which is the 
column density for which the typical optical depth for the photon to escape is unity). 
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Some photons are also emitted by helium recombinations. He i makes little difference to 
the total budget of ionizing photons, because for each absorption by Hei , one photon that 
ionizes hydrogen is produced by the subsequent recombination. Because we did not include 
absorption by Hei in our previous estimate of the mean free path, we do not include the 
reemission either. 


Finally, photons with energy above 4 Rydbergs are mostly used to photoionize Hen , 
and each He ii recombination can produce between one and three photons capable of ionizing 
hydrogen. About 25% of the recombinations that are not direct to the ground state go to 
n = 2, producing a Balmer continuum photon, which can ionize hydrogen. In addition, all 
the recombinations produce either a Lyman series photon, or two continuum photons when 
the last transition is from 2 s to Is; in ~ 8 % of the recombinations, two continuum photons 
that can both ionize hydrogen are produced (see Osterbrock 1989, Tables 2.8, 4.11, and 4.12). 
Hence, on average 1.33 ionizing photons are produced for each Hen that is photoionized. For 
a source spectrum oc with rj = 1.5, the fraction of photons used to ionize Hen is 

1/8, so Hen reemission increases the emissivity by a factor 1 + 0.33/8 = 1.04. This factor 
can be reduced during the Hen reionization, when the rate of recombinations is less than 
the rate of photoionizations. 


Based on this discussion of the effect of reemission from absorbers, we reduce our upper 
limit in ( 6 ) for the total emissivity by a factor 1.14, to obtain the upper limit for the direct 
emissivity from sources at z = 4 (rounded up to one digit given the various uncertainties 
discussed). 
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3. Comparison with the emissivity from Lyman break galaxies and QSOs 

3.1. QSOs 

Fan et ah (2001) measured the following QSO luminosity function at high redshift 
{z > 3.5) from SDSS; 

iI)q{L) dL = 7.8 X 10"® Mpc"® f 1 q-o.47 (^-3)^^ ^ 

V-h-26/ -h_26 

where L_ 2 q is the luminosity of a QSO with absolute AB magnitude M 1450 = —26 at A = 
1450 A, and 7 = 2.58. Using a power-law spectral index rjQ = 1.5 for the QSO spectrum at 
A < 1450 A (Telfer et al. 2002 hnd rjq = —1.57 at A < 1200 Afor radio-quiet QSOs, and the 
spectral slope is shallower at longer wavelengths; e.g., Schmidt, Schneider, & Gunn 1995), the 



total number of ionizing photons (with A < 912 A) emitted is L- 2 & = 5.48 x 10^®photons/s. 
The luminosity function in equation (8) was observed down to L ~ 0.5L_26 by Fan et 
ah (2001); the total emissivity from QSOs clearly depends on the luminosity at which the 
slope of 'ipqiL) becomes shallower than L~‘^. If we assume the form 'ipq^L) oc [{L/Lq^)'^~^ + 
(L/Lq*)'^]“\ which produces a similar slope of t/’(F) at low luminosities {L < L*) as observed 
at lower redshift (e.g., Boyle et ah 1988, 2000), then the comoving emissivity is 
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Expressing this as the number of photons emitted per Hubble time and per atom, at 
we hud. 


^qso 

H{z = A)nat 



(9) 
= 4, 

( 10 ) 


A lower limit to the emissivity from QSOs is obtained if L* ~ 0.5L_26, ^qsol{.H^at) ^ 1- 
In order not to exceed the upper limit obtained in the last section, e/[H{z = ^)nat] < 7, 
we must require Lq* > 0.015L_26. We note, however, that low-luminosity QSOs may be 
intrinsically self-absorbed at the Lyman limit, reducing their contribution to the ionizing 
emissivity (see Alam & Miralda-Escude 2002). 


The number of faint AGNs has been constrained by Steidel et ah (2002), who hud 
that their abundance at z = 3 and at an apparent AB magnitude R ~ 23.9 (or luminosities 
L = 0.018L_26) is about 3% of the Lyman break galaxies abundance at a magnitude R = 24.6 
at 2 : = 3. Using the Lyman break galaxy luminosity function given in the next subsection, 
we hud that this implies that the turnover of the QSO luminosity function should occur at 
Fq* = 0.046L_26, assuming again the form tjjq oc [(L/Lq*)'''”^ -|- (L/Lq*)'’']”^. If the value 
of Lq* were the same at = 4 as at 2 : = 3 (but with the QSO abundance having declined 
as 10“°-"^’^^ at all luminosities, as in eq. [8]), QSOs would contribute about half of the total 
emissivity required at = 4 from equation (7). There are also limits from the Hubble Deep 
Field (Conti et al. 1999), but these probe QSOs at much lower luminosities and do not 
provide a stronger constraint on Lg*. In summary, the available observations are consistent 
with a contribution to the ionizing emissivity from QSOs that does not exceed the upper 
limit in (7). 
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3.2. Lyman break galcixies 

The galaxy ultraviolet luminosity function at high redshift has been measured by Steidel 
et ah (1999). They £t the luminosity function to a Schechter function, 

( 11 ) 

After converting to the flat cosmological model with = 0.3, h = 0.65, the parameters they 
hnd are 0* = 8.6 x 10“"^Mpc“^, a = 1.6, and L* = 1.4 x 10^® erg s“^ Hz“^ at A = 1700 A. To 
convert this to the ionizing luminosity, we use the models of Bruzual & Chariot (1993), for 
which the luminosity per unit frequency decreases by a factor 10 from A = 1700 A to 912 A, 
for an age of 10® years (these models assume a Salpeter mass function up to a maximum mass 
of 125M0 and solar metallicity). Using also the approximation of L,y oc at A < 912 A, 
we hnd the ionizing luminosity is L* = 1.06 x 10®"^ photons/ s. 

The form of the galaxy luminosity function (11) is highly uncertain at L < L*, because 
the observations of Steidel et al. (1999) are not deep enough to reach lower luminosities, and 
galaxies down to 10“® ®L* are only detected in the Hubble Deep Field (e.g., Madau et al. 
1996). The Hubble Deep Field counts at z = 4 are significantly below the prediction from the 
model in equation (11), but this may be due to galaxy clustering (Steidel et al. 1999; see their 
Fig. 8). If we include only galaxies with luminosities larger than (L*, 10 — 0.4L*, 10“® ®L*), 
the total emissivity inferred from (11) is 

^ = (3-4, 8.5,13.8)/e.c [L > (1,10-®•^ 10-o-«)L,] . (12) 

Steidel, Pettini, & Adelberger (2001) found that Lyman break galaxies may have escape 
fractions of ionizing photons as large as 50%, from a possible detection of Lyman contin¬ 
uum hux after coadding spectra of several galaxies. If this result were true, the luminous 
galaxies with L > 10“® ®L* would already produce the maximum emissivity in equation (7). 
Alternatively, fesc may be much lower (see Giallongo et al. 2002; note also that the value 
fesc — 0.5 measured by Steidel et al. 2001 was for the subset of bluest quartile of the Lyman 
break galaxies, which may not be representative), and lower luminosity objects may be re¬ 
sponsible for the bulk of the emissivity. More accurate determinations of the escape fraction, 
the abundance of Lyman limit systems, and the inferred photoionization rate from the Lya 
forest transmitted hux will be needed to resolve this question. 
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4. Constraints on the emissivity evolution to end reionization at z = 6 

The highest redshift QSO known at the present time (Becker et ah 2001) shows a 
complete Gunn-Peterson trough at 2: ~ 6, but the spectra of all the sources at lower redshift 
have detectable transmitted flux blueward of Lya , indicating that most of the volume in the 
intergalactic medium had to be reionized hy z = 6. We now address the contraint implied 
by this fact on the evolution of the ionizing emissivity, in view of its value at z = 4 (eq. 7). 


4.1. The reionization equation 

Reionization proceeds by the creation of ionized regions around individual sources 
(Arons & Wingert 1972). The mean free path of ionizing photons through the neutral 
medium is typically very short compared to the distance between sources, so reionization 
proceeds by the growth of the ionized regions and can be modeled in terms of the fraction 
of the volume that is ionized at time t, y{t), assuming that the remaining fraction 1 — y is 
completely neutral. Moreover, because the emitted photons are used to ionize atoms almost 
immediately after their emission, it follows that the change in the fraction of atoms that 
are ionized in the universe is equal to the ionizing photon emission rate per atom minus the 
average recombination rate of atoms. If ea{t) = €i{t)/nat is the comoving ionizing photon 
emissivity per atomic nucleus, R{t) is the average recombination rate per atom in ionized 
regions (which includes the effect of an effective clumping factor of the ionized gas), and T) 
is the fraction of baryons that actually need to be ionized in the ionized fraction y of the 
universe (with the remaining 1 — T) fraction being the baryons that are in dense systems, 
such as stars or dense gas in galaxies, that do not need to be ionized by external sources), the 
equation for the evolution of y is (Madau, Haardt, & Rees 1999; Miralda-Escude, Haehnelt, 
& Rees 2000) 

d{yFi) 


dt 


= e„ 


Ry . 


(13) 


This equation (apart from the factor T)) is analogous to the one obtained for the problem 
of an expanding cosmological Hii region around an ionizing source, which was discussed by 
Shapiro & Giroux (1987) and Donahue & Shull (1987), where y behaves as the volume of 
the Hii region produced by a source of luminosity proportional to eait). It is also identical 
to the equation of radiative transfer, where yFi plays the role of the intensity, and RF~^ 
and Ca play the role of the absorption and emission coefficients, respectively. The solution is 
found through the substitutions dr = RF~^ dt, and w = yFi e^\ 


y{t) = F-^ / , 


(14) 
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where r(t) = f* dt"F~^R{t"), and U is the initial time when the first sonrces turn on. 
Assuming that any clumping factor is constant, and that — 1 [e.g., for the cosmological 

constant model with a present matter density hlmo = 0-3, we have Vtm{z = 4) = 0.982], the 
recombination rate varies as R{t) oc Choosing = t{z = 4) = as a fiducial 

time, and R 4 = R{t 4 ), we have: 

Vit) = d{t'/U) exp R 4 F-H 4 • (15) 

*4 

4.2. The recombination rate 

Before presenting solutions of the reionization equation, we discuss the value that we 
adopt for the recombination rate, R. The average recombination rate is affected by the 
clumping factor of the ionized gas, which we dehne as: Ci =<nl>j / <ne>‘^j. The subscript 
J is used here to denote that the average is to be taken only over those regions of space where 
matter is predominantly photoionized by photons that have escaped to the intergalactic 
medium. This restriction is necessary because if we did not impose it, the clumping factor 
would be dominated by regions of dense gas that are locally ionized by sources that were not 
included in the emissivity in equation (6), precisely because the radiation from these sources 
is locally absorbed. In fact, in the absence of any restriction on the regions over which one 
must average to obtain the clumping factor of ionized matter, then the clumping factor would 
obviously be very large and dominated by the densest stars. Once stars are eliminated from 
consideration, there are still H ii regions that form in dense clouds in the interstellar medium 
of galaxies around massive young stars, in which all the emitted photons are immediately 
absorbed locally. This emission is irrelevant for the reionization of the intergalactic medium, 
and it is best not to include in equation (6), and at the same time not to include the dense 
gas ionized locally in the clumping factor. This dehnition of the clumping factor inevitably 
involves an arbitrary choice in defining when a photon is considered to be escaped, and 
therefore when a region is predominantly ionized by escaped photons rather than local ones. 
In practice, a photon can be dehned to have escaped once it reaches an intergalactic region 
with gas density below some critical value. The calculation of the clumping factor can only 
be obtained from full radiative transfer simulations of reionization, but its value can differ 
among studies by large factors depending on the dehnition that is adopted (see Gnedin 
2000 and references therein). If the restriction that escaped photons dominate the radiation 
intensity is not imposed, then very large values of the clumping factor can be obtained, which 
could only be used consistently by including the locally absorbed photons in the emissivity. 

Here, we use the simple model of Miralda-Escude et ah (2000) for the clumping factor. 
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which uses the approximation that the gas is ionized only up to an overdensity Aj = Pi/p, 
and all the gas at higher densities is neutral, or otherwise is ionized locally and therefore 
does not contribute to the clumping factor. The overdensity Aj is related to the mean free 
path Aj that a photon can traverse before being absorbed. During reionization, the mean 
free path A* should be of order the size of the H ii regions, or of order the typical separation 
between neighboring sources. For A* = 1000 km s~^, or ~ lOMpc aX, z = Qin comoving units, 
the critical overdensity is Aj ~ 5, and the clumping factor is very close to unity (see Figs. 
2d and 6 in Miralda-Escude et ah 2000). 

It has been pointed out by Haiman, Abel, & Madau (2001) that before reionization 
is complete, a large number of low-mass halos should be present and should contain dense 
gas accreted before it was heated by reionization. The effect of these low-mass halos would 
be to decrease the mean free path at a hxed overdensity Aj, and to widen the distribution 
of the overdensity. To estimate the plausible change in the clumping factor due to these 
low-mass halos, we can assume that the structure of the atomic intergalactic medium before 
reionization is self-similar to that of the ionized intergalactic medium, but with the scales 
reduced in proportion to the Jeans scale Aj, so that the clumping factor should depend only 
on the ratio Aj/Aj, and the amplitude of the primordial density fluctuations at the Jeans 
scale. If the temperature increases by a factor 100 during reionization (from ~ 100 K to 
~ 10^ K), the Jeans length increases by a factor 10. The amplitude of primordial density 
fluctuations in Cold Dark Matter with VLm = 0.3 and h = 0.65 increases by a factor ~ 2 
between the scales of 0.05 to 0.5 Mpc (roughly corresponding to the Jeans scales at these two 
temperatures). Therefore, the clumping factor of a medium in which all the low-mass halos 
formed from scales intermediate between the Jeans scales of the atomic and ionized medium 
have survived, at z = 7 and \i = lO^kms”^, should be the same as the clumping factor 
of a medium with no such surviving halos, evaluated at z = 3 and Aj = lO^kms”^. From 
Figures 2b and 6 of Miralda-Escude et al. (2000), we hnd this clumping factor to be 7. We 
note that this is an upper limit to the true clumping factor, because the gas in the low-mass 
halos that might increase the clumping factor above the value of ~ 1 found previously is 
likely to evaporate on a time short compared to the Hubble time, once these halos have been 
reached by an ionization front (Shapiro & Raga 2001). The gas may also be expelled by a 
small amount of internal star formation even before the halo is reached by external ionizing 
radiation. 

Finally, a discussion is needed to decide if the case A or case B recombination coefficient 
should be used. Ionizing photons of frequency u will typically be absorbed in regions of 
column density Nfji ~ 1-6 x 10^^ where uhi is the Lyman limit frequency. The 

sources typically have power-law spectra and most of the emitted photons have frequencies 
substantially above uhi, especially if the spectra are hardened by internal absorption. Most 
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of the photons emitted by direct recombinations to the ground state, which have frequencies 
very close to the Lyman limit, will therefore be produced in high column density systems 
and will be reabsorbed within the same system (in addition, among the small fraction of 
photons from direct recombinations to the ground state that escape the system in which 
they have been produced, some will be redshifted below the Lyman limit and will not be 
absorbed again). Hence, these photons do not contribute to increase the cosmic ionizing 
background intensity for the most part. An alternative way to understand this argument is 
that, even though in principle we are in a case B situation because the recombination photons 
are reabsorbed, in practice the inclusion of the photons from recombinations increases the 
clumping factor by the ratio of the case A to case B recombination coefficients because these 
photons increase the ionization in dense regions, therefore yielding a recombination rate 
equivalent to using case A. Therefore, we use the case A recombination coefficient. 


4.3. Results for the reionization history 

We now compute y{t) assuming models where the emissivity varies as Ca oc (1 + z)~^ oc 

^ 2 / 3/3 

z > and has the value ea/H 4 = 7 at = 4, i.e., the maximum allowed value found 
previously (eq. 7). We choose the recombination coefficient a = 3.15 x 10“^^cm^s“^, valid 
for case A and T = 1.5 x 10^ K, and electron density Ue = 6.95 x 10“® cm“^ at x = 6, valid 
for = 0.022 and if helium is only once ionized. This yields a recombination rate per 
Hubble time R{z = 6 )/H{z = 6 ) = 1.02. In some cases, we increase this recombination rate 
by a clumping factor Ci. We £x the fraction of baryons in the ionized regions that actually 
need to be ionized to T) = 0.9 (independent of redshift), in reasonable agreement with the 
fraction of baryons that seem to be present in damped Lya systems at x ~ 4 (Lanzetta et al. 
1991; Storrie-Lombardi et al. 1996), and the fraction of baryons above the critical overdensity 
Aj ~ 5 at = 6 in the model discussed above for the clumping factor (Miralda-Escude et 
al. 2000, Fig. 2d). 

We solve equation (15), starting at y = 0 at some initial redshift x* when the emissivity 
is assumed to turn on suddenly. The parameters of the model that we vary are Zi, (3, and 
Cl, and we require all models to reach y = I (i.e., the end of reionization) at x = 6. 

The results of five representative models are shown in Figure 1. When Ci = 1, assuming 
the emissivity to remain constant satisfies the requirement of completing reionization by 
z = 6 when the emissivity is turned on at Zi = 8.7. With only a slight decline of the 
emissivity, the initial redshift needs to be pushed to much higher values. For (3 = 0.63, 
which corresponds to a decline of the comoving emissivity by a factor of only ~ 1.5 between 
z = A and z = 9, the sources need to start at very high redshift to reach y = 1 at z = 6 . 




Fig. 1.— Fraction of the volume of the universe Filed by ionized regions as a function of 
redshift. The models assume the emissivity is equal to ea/H{t) = 7[(1 + z)/h)]~^^ at z < Zi, 
where Zi is the initial redshift where the emission is turned on. The clumping factor of the 
ionized gas is Ci. The models are constrained to reach a fraction of unity at z = 6. The 
electron scattering optical depth to CMB photons is indicated for each model. 
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Naturally, changes in the emissivity at redshifts higher than ~ 15 are irrelevant for the 
requirement that reionization ends at z = 6, because most of the atoms ionized at this high 
redshift will recombine again. We therefore conclude that, for the case of no clumping, the 
ionizing emissivity at 6 < z < 9 cannot he much lower than its maximum value allowed at 
z = 4. In other words, the sources responsible for the emission of ionizing radiation cannot 
decrease their overall comoving emissivity by more than a factor ~ 1.5 between z = 4 and 
= 9. 

Allowing now for a clumping factor, an increase to Ci = 2.1 already requires the model 
with constant emissivity {(3 = 0) to maintain the emissivity up to a very high initial redshift. 
A model where the same amount of radiative energy is emitted at every Hubble time [i.e., 
f3 = —1.5, since H{z) oc (1 + 2:)^'®], starting at very high 2:^, results in complete reionization 
at z = 6 when Ci = 5.8. Therefore, for a clumping factor as large as ~ 7, as suggested 
previously for a model where the ionized gas in all low-mass halos is retained, the emissivity 
would need to substantially increase with redshift. 


4.4. Comparison to the observed evolution at 2: < 4 

We have shown that by combining the measurement of the ionizing emissivity e* at 
z = 4 with the requirement that the reionization must end by = 6, we can infer that 
the emissivity does not decrease with redshift by more than a factor ~ 1.5 from z = 4 to 
z = 9. We can now ask how this evolutionary constraint compares to the evolution at z < 4 
that is implied by observations of the Lya forest decrement and the mean free path between 
Lyman limit systems. From McDonald & Miralda-Escude (2001), the photoionization rate 
(or proper intensity of the background) declines roughly as (1 -|- z)~^'^ from 2: = 2.5 to 
z = 4. The abundance of Lyman limit systems per unit redshift is proportional to (1 -f z)^'^ 
(Storrie-Lombardi et ah 1994), implying that the proper mean free path varies as (1 -|- z)~'^. 
Hence, the comoving emissivity evolves as (1 -|- z)~^'^, essentially constant within the likely 
errors of measurement for both the proper intensity and mean free path evolution. Therefore, 
the results of §4.3 show that if the clumping factor is low, this constant emissivity needs to 
roughly maintained at the same value up to 2 ~ 9, while if the clumping factor is much larger 
than one, the emissivity needs to modify its evolution and increase strongly with redshift at 
z > 4. 
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4.5. The Thomson optical depth of the Cosmic Microwave Background 


The optical depth of Cosmic Microwave Background photons to electron scattering, Tes, 
is given by 

r^i c dt 

Tes ^ j dz + zfaes y{z) Fi , 


(16) 


where (7^.^ is the Thomson cross section for electron scattering. The values of Tes for the five 
models in Figure 1 are indicated in the hgure. We have assumed Fi = 0.9 at all redshifts; 
the 10% of the matter not included is supposed to be either gas that remains atomic or 
molecular in self-shielded regions, or stars. 


These results show that the minimum allowed value of Tes, which occurs for models where 
the emissivity turns on at a redshift not much larger than 6, is ~ 0.045. The maximum values 
of the optical depth depend on the emissivity at high redshift. If the emissivity does not 
increase with redshift, then the ionized fraction of the volume at high redshift needs to be 
small due to the fast rate of recombinations (i.e., for constant emissivity, ?/ oc (1 -|- z)~^, and 
the integral in eq. [16] converges rapidly at high z). For a model with a constant emissivity 
up to Zi = 20 and clumping factor equal to one, the optical depth can be increased only up 
to Tes = 0.09 (in this model, reionization ends at z = 7.3, but the reionization end could 
always be delayed by increasing the clumping factor as the fraction y reaches a value close 
to unity). 

An emissivity increasing with redshift can obviously produce larger optical depths, and 
as mentioned before the value of the emissivity at z >15 does not affect the epoch at which 
reionization ends, so in principle the optical depth could be arbitrarily high if enough ionizing 
sources were present at very high redshift. However, in the flat Cold Dark Matter model 
with Qrn = 0.3, which is strongly supported by all observations of large-scale structure, 
the hrst stars form only at ~ 20 from 3 — a peaks collapsing into the hrst low-mass halos 
(e.g., Couchman & Rees 1986; Bromm, Coppi, & Larson 1999; Abel, Bryan, & Norman 2000; 
Ricotti, Gnedin, & Shull 2002; Venkatessan, Tumlinson, & Shull 2003), and so the emissivity 
cannot continue to increase up to a very high redshift. At the same time, observations 
indicate a cosmic star formation rate in galaxies that starts to decline with redshift at z >2 
(Madau, Pozzetti, & Dickinson 1998), although the ionizing emissivity might increase with 
redshift even when the star formation rate decreases because of an increasing escape fraction 
for ionizing photons. 


As this paper was being refereed, the results of the WMAP mission were announced. 
Kogut et ah (2003) give a result Tg = 0.16±0.04 from a model-independent analysis, although 
the result and the error can change when htted to different CDM models (Spergel et al. 2003). 
The value of ~ 0.08 implied by a constant emissivity with clumping factor of one up to very 
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high redshift would represent a 2 — a deviation from the value of Kogut et ah (2003). A 
value as high as = 0.16 would clearly imply a substantial increase of the emissivity up to 
~ 20, when only a small fraction of the mass has collapsed into halos that can form stars 
in the CDM model. 


5. Conclusions 

The observed mean transmitted flux of the Lyo forest, when combined with indepen¬ 
dently measured values of the mean baryon density, the amplitude of the power spectrum, 
and the gas temperature, yields a measurement of the intensity of the ionizing background. 
The abundance of Lyman limit systems determines the mean free path of ionizing photons. 
These two quantities together determine the ionizing emissivity. We have argued in this 
paper that the emissivity measured in this way at 2: = 4 is at most 7 ionizing photons per 
atom per Hubble time. We have then shown that the emissivity over the redshift range 
6 < 2: < 9 cannot be much lower than the value at = 4 in order that the universe can 
be reionized hy z = 6, if the clumping factor of ionized gas during the reionization epoch is 
close to unity. If the clumping factor is larger than unity, then the ionizing emissivity must 
increase with redshift between z = 4 and z ^ 9. 

The electron scattering optical depth must be at least 0.045, given the presence of ionized 
intergalactic gas at z < 6, and a reasonable minimum redshift range over which sources ionize 
the universe at z > 6. If the emissivity of ionizing photons does not increase with redshift 
at z > 4, then the optical depth must be less than 0.09. Measurements of the optical depth 
from the CMB spectrum of temperature and polarization fluctuations therefore provide a 
powerful diagnostic to decide if the peak of the comoving ionizing emissivity occurred at 
z 4 (with a broad plateau extending up to ~ 9 so that the end of reionization can be 
reached by = 6), or at a much earlier epoch. The recent WMAP result has an error that is 
still too large to reach a dehnite conclusion, but if a value Tg — 0.16 is conhrmed it implies 
that the comoving emissivity continues to rise up to ~ 20. 


This work was supported in part by grant NSF-0098515. 
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